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ABSTRACT 

Recent estimates of cosmological parameters derived from Cosmic Microwave Back- 
ground (CMB) anisotropics are based on the assumption that we know the precise 
amount of energy density in relativistic particles in the universe, uj r ei, at all times. 
There are, however, many possible mechanisms that can undermine this assumption. 
In this paper we investigate the effect that removing this assumption has on the de- 
termination of the various cosmological parameters. We obtain fairly general bounds 
on the redshift of equality, z eq = ui m /oj re i = 3100^4qq. We show that u> re i is nearly 
degenerate with the amount of energy in matter, w m , and that its inclusion in CMB 
parameter estimation also affects the present constraints on other parameters such as 
the curvature or the scalar spectral index of primordial fluctuations. This degeneracy 
has the effect of limiting the precision of parameter estimation from the MAP satellite, 
but it can be broken by measurements on smaller scales such as those provided by the 
Planck satellite mission. 
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1 INTRODUCTION 

Our knowledge of the cosmological parameters has increased 
dramatically with the release of recent cosmic microwave 
background (CMB) observations (Netterfield et al. 2001; Lee 
et al. 2001; Halverson et al. 2001). Recent analyses of these 
datasets (de Bernardis et al. 2001; Pryke et al. 2001; Wang, 
Tegmark & Zaldarriaga 2001) have reported strong new con- 
straints on various parameters including the curvature of 
the universe and the amount of baryonic and dark matter. 
The precise determination from the CMB of other parame- 
ters such as the cosmological constant or the spectral index 
of primordial fluctuations can be limited by various degen- 
eracies, and such degeneracies are best lifted by combining 
CMB data with either supernova (SN) data (Garnavich et al. 
1998; Perlmutter et al. 1997) or large scale structure (LSS) 
surveys, such as PSCz (Saunders et al. 2000), 2dF (Percival 
et al. 2001) or Lyman-a (Croft et al. 2001) data. At present, 
the values obtained from CMB measurements under the as- 
sumption of purely adiabatic fluctuations are consistent with 
the generic predictions of the inflationary scenario, n s ~ 1 
and Q,tot = 1 (Linde 1990), and with the standard big-bang 
nucleosynthesis bound, Q.bh 2 = 0.020 ± 0.002 (Buries et al. 
2001; Esposito et al. 2001). All of these observations con- 



verge towards a consistent picture of our universe, provid- 
ing strong support for the inflationary scenario as the mech- 
anism which generated the initial conditions for structure 
formation. 

The derivation of the cosmological parameters from 
CMB is, however, an indirect measurement, relying on the 
assumptions of a theoretical scenario. For this reason, recent 
efforts have been made to study the effects of the removal 
of some of these assumptions. For example, a scale-invariant 
background of gravity waves, generally expected to be small, 
has been included in the analysis of Kinney, Melchiorri & Ri- 
otto (2000), Wang et al. (2001) and Efstathiou et al. (2001), 
with important consequences for parameter estimation. A 
scale-dependence of the spectral index has been included in 
the analysis of Griffiths, Silk & Zaroubi (2001), Santos et al. 
(2001) and Hannestad et al. (2001). Furthermore, in Bucher, 
Moodley & Turok (2000), Trotta, Riazuelo & Durrer (2001) 
and Amendola et al. (2001), the effects of including isocur- 
vature modes, which naturally arise in the most general in- 
flationary scenarios, have been studied, with the finding that 
the inclusion of these modes can significantly alter the CMB 
result. Even more drastic alterations have been proposed in 
Bouchet et al. (2001) and Durrer, Kunz & Melchiorri (2001). 

All the above modifications primarily affect the con- 
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straints on the curvature, on the physical baryon density 
parameter, u)b = flbh 2 , and the scalar spectral index n s . 

In this paper we study another possible modification 
to the standard scenario, namely variations in the param- 
eter uj r ei which describes the energy density of relativistic 
particles at times near decoupling, T ~ 0.1 eV. CMB data 
analysis with variations in this parameter has been recently 
undertaken by many authors (Hannestad 2000; Esposito et 
al. 2001; Orito et al. 2001; Hansen et al. 2001; Kneller et 
al. 2001; Hannestad 2001; Zentner & Walker 2001), giving 
rather crude upper bounds, significantly improved only by 
including priors on the age of the universe or by includ- 
ing supernovae (SN) or large scale structure (LSS) data. It 
is worth emphasizing that there is little difference in the 
bounds obtained on N e ff, the effective number of relativis- 
tic species, between old and recent CMB data because of the 
degeneracy which we will describe in detail below. We focus 
here on the effects that the inclusion of this parameter, u> r ei, 
has on the constraints of the remaining parameters in the 
context of purely adiabatic models. 

As we will show below (and as observed previously, see 
e.g. Hu et al. (1999)) there is a strong degeneracy between 
ui re i and io m . This is important, because an accurate de- 
termination of uj m from CMB observations (and of fl m by 
including the Hubble Space Telescope result h — 0.72±0.08) 
can be useful for a large number of reasons. First of all, de- 
termining <jj c dm = Mm — <jJb can shed new light on the nature 
of dark matter. The thermally averaged product of cross- 
section and thermal velocity of the dark matter candidate 
is related to u> m , and this relation can be used to analyze 
the implications for the mass spectra in versions of the Su- 
persymmetric Standard Model (see e.g. Barger & Kao 2001, 
Djouadi, Drees & Kneur 2001, Ellis, Nanopoulos & Olive 
2001). The value of Q m can be determined in an indepen- 
dent way from the mass-to-light ratios of clusters (Turner 
2001), and the present value is 0.1 < Q m < 0.2 (Carlberg et 
al. 1997; Bahcall et al. 2000). Furthermore, a precise mea- 
surement of Q m will be a key input for determining the red- 
shift evolution of the equation of state parameter w(z) and 
thus discriminating between different quintessential scenar- 
ios (see e.g. Weller and Albrecht 2001). 

This paper is structured as follows. In the next section, 
we briefly review various physical mechanisms that can lead 
to a change in io re i with respect to the standard value. In 
section 3, we illustrate how the CMB angular power spec- 
trum depends on this parameter and identify possible de- 
generacies with other parameters. In section 4, we present 
a likelihood analysis from the most recent CMB data and 
show which of the present constraints on the various param- 
eters are affected by variations in uj re i. Section 5 forecasts 
the precision in the estimation of cosmological parameters 
for the future space missions MAP and Planck. Finally, in 
section 6, we discuss our results and present our conclusions. 



2 EFFECTIVE NUMBER OF RELATIVISTIC 
SPECIES 

In the standard model cu re i includes photons and neutrinos, 
and it can be expressed as 

LU re l = U 1 + Neff ■ LOu (1) 



where oj 7 is the energy density in photons and uj v is the 
energy density in one active neutrino. Measuring u re i thus 
gives a direct observation on the effective number of neutri- 
nos, N e ff. Naturally there are only 3 active neutrinos, and 
Neff is simply a convenient parametrization for the extra 
possible relativistic degrees of freedom 

N eff = 3 + ANcmb . (2) 

Thus u) re i includes energy density from all the relativis- 
tic particles: photons, neutrinos, and additional hypothet- 
ical relativistic particles such as a light majoron or a ster- 
ile neutrino. Such hypothetical particles are strongly con- 
strained from standard big bang nucleosynthesis (BBN), 
where the allowed extra relativistic degrees of freedom typ- 
ically are expressed through the effective number of neutri- 
nos, N e ff = 3 + ANbbn- BBN bounds are typically about 
ANbbn < 0.2 - 1.0 (Buries et al. 1999; Lisi, Sarkar & Vil- 
lante 1999). 

One should, however, be careful when comparing the 
effective number of neutrino degrees of freedom at the times 
of BBN and CMBR, since they may be related by differ- 
ent physics (Hansen et al. 2001). This is because the energy 
density in relativistic species may change from the time of 
BBN (T ~ MeV) to the time of last rescattering (T ~ eV). 
For instance, if one of the active neutrinos has a mass in 
the range eV < m < MeV and decays into sterile par- 
ticles such as other neutrinos, majorons etc. with lifetime 
t(BBN) < r < i(CMBR), then the effective number of neu- 
trinos at CMBR would be substantially different from the 
number at BBN (White, Gelmini & Silk 1995). Such mas- 
sive active neutrinos, however, do not look too natural any 
longer in view of the recent experimental results on neutrino 
oscillations (Fogli et al. 2001; Gonzalez-Garcia et al. 2001), 
showing that all active neutrinos are likely to have masses 
smaller than an eV. One could instead consider sterile neu- 
trinos mixed with active ones which could be produced in 
the early universe by scattering, and subsequently decay. 
The mixing angle must then be large enough to thermal- 
ize the sterile neutrinos (Langacker 1989), and this can be 
expressed through the sterile to active neutrino number den- 
sity ratio n s /n„ « 4-10 4 sin 2 29 (m/keV)(10.75/3*) 3/2 (Dol- 
gov & Hansen 2001), where 9 is the mixing angle, and g* 
counts the relativistic degrees of freedom. With n s /n„ of or- 
der unity we use the decay time, r ~ 10 20 (keV/m) 5 / sin 2 29 
sec, and one finds, r w 10 17 (keV/m) 4 yr, which is much 
longer than the age of the universe for m ~ keV, so they 
would certainly not have decayed at t(CMBR). A sterile 
neutrino with mass of a few MeV would seem to have the 
right decay time, r ~ 10 5 yr, but this is excluded by stan- 
dard BBN considerations (Kolb et al. 1991; Dolgov, Hansen 
& Semikoz 1998). 

Even though the simplest models predict that the rela- 
tivistic degrees of freedom are the same at BBN and CMB 
times, one could imagine more inventive models such as 
quintessence (Albrect & Skordis 2000; Skordis & Albrect 
2001) which effectively could change AN between BBN and 
CMB (Bean, Hansen & Melchiorri 2001). Naturally AN can 
be both positive and negative. For BBN, AiV can be nega- 
tive if the electron neutrinos have a non-zero chemical po- 
tential (Kang & Steigman 1992; Kneller et al. 2001), or more 
generally with a non-equilibrium electron neutrino distribu- 
tion function (see e.g. Hansen & Villante 2000). To give an 
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explicit (but highly exotic) example of a different number 
of relativistic degrees of freedom between BBN and CMB, 
one could consider the following scenario. Imagine another 
2 sterile neutrinos, one of which is essentially massless and 
has a mixing angle with any of the active neutrinos just big 
enough to bring it into equilibrium in the early universe, and 
one with a mass of rn Vs = 3 MeV and decay time t„ 3 =0.1 
sec, in the decay channel v s — *• + cj>, with <f> a light scalar. 
The resulting non-equilibrium electron neutrinos happen to 
exactly cancel the effect of the massless sterile state, and 
hence we have ANbbn = 0. However, for CMB the picture 
is much simpler, and we have just the stable sterile state 
and the majoron, hence ANomb = 1.57. For CMB, one can 
imagine a negative AN from decaying particles, where the 
decay products are photons or electron/positrons which es- 
sentially increases the photon temperature relative to the 
neutrino temperature (Kaplinghat & Turner 2001). Such a 
scenario naturally also dilutes the baryon density, and the 
agreement on uit from BBN and CMB gives a bound on 
how negative ANcmb can be. Considering all these possi- 
bilities, we will therefore not make the usual assumption, 
ANbbn = ANcmb, but instead consider ANcmb as a 
completely free parameter in the following analysis. 

The standard model value for N e f j with 3 active neutri- 
nos is 3.044. This small correction arises from the combina- 
tion of two effects arising around the temperature T ~ MeV. 
These effects are the finite temperature QED correction to 
the energy density of the electromagnetic plasma (Heck- 
ler 1994), which gives AN = 0.01 (Lopez & Turner 1999; 
Lopez et al. 1999). If there are more relativistic species 
than active neutrinos, then this effect will be correspond- 
ingly higher (Steigman 2001). The other effect comes from 
neutrinos sharing in the energy density of the annihilating 
electrons (Dicus et al. 1982), which gives AN = 0.034 (Dol- 
gov, Hansen & Semikoz 1997, 1999; Esposito et al. 2000). 
Thus one finds iV e // = 3.044. It still remains to accurately 
calculate these two effects simultaneously. 



3 CMB THEORY AND DEGENERACIES 

The structure of the Ce spectrum depends essentially on 4 
cosmological parameters 



, uirei and TZ , 



(3) 



the physical baryonic density u;, = Q,bh 2 , the energy den- 
sity in matter w m = (f2 C£ j m + Qb)h 2 , the energy density in 
radiation u re i and the 'shift' parameter TZ = £ re f/£, which 
gives the position of the acoustic peaks with respect to a flat, 
SIa = reference model. Here h denotes the Hubble param- 
eter today, Ho = WOh km s _1 Mpc -1 , and Qa is the density 
parameter due to a cosmological constant, Oa = A/3Hq. In 
previous analyses (Efstathiou & Bond 1999; Melchiorri & 
Griffiths 2000 and references therein), the parameter ui re i 
has been kept fixed to the standard value, while here we 
will allow it to vary. It is therefore convenient to write 
to rel =4.13- 10~ 5 (1 + 0. 135 -AiV) (taking Toms = 2.726 K), 
where A./V is the excess number of relativistic species with 
respect to the standard model, AT e // = 3 + A./V. The shift 
parameter TZ depends on Q m = Q, c dm on the curvature 

flk and on f2 re ; = uj re i/h 2 through 
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Fig. 1 — Top panel: CMB degeneracies between cosmological 
models. Keeping z e q,iOb an <i TZ fixed while varying AN produces 
nearly degenerate power spectra. The reference model (black, 
solid) has AN = 0, Qtot = 1.00, n s = 1.00; the nearly degen- 
erate model (blue, dotted) has AN = 10, Q to t = 1-05, n s = 1.00. 
The position of the peaks is perfectly matched, only the relative 
height between the first and the other acoustic peaks is some- 
what different in this extreme example. The degeneracy can be 
further improved, at least up to the third peak, by raising the 
spectral index to n s = 1.08 (red, dashed). Bottom panel: the 
matter power spectra of the models plotted in the top panel to- 
gether with the observed decorrelated power spectrum from the 
PSCz survey (Hamilton and Tegmark 2000). The geometrical de- 
generacy is now lifted. 




[fUi(l + zf + fim(l + Z) A + O fe (l + zy + 



(4) 



(5) 



i-l/2 



The function x(v) depends on the curvature of the uni- 
verse and is y, sin(y) or sinh(j/) for flat, closed or open mod- 
els, respectively. Eq. (^) generalizes the expression for TZ 
given in Melchiorri & Griffiths (2001) to the case of non- 
constant Q re i. 

By fixing the 4 parameters given in Q), or equivalently 
the set u)b, the redshift of equality z eq = uj m /u re i, AN and 
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Fig. 2. — H as a function of AN with Ha = 0.7 and O m = 0.3. 
The position of the peaks is only weakly affected by AN. 

1Z, one obtains a perfect degeneracy for the CMB anisotropy 
power spectra on degree and sub-degree angular scales. On 
larger angular scales, the degeneracy is broken by the late 
Integrated Sachs- Wolfe (ISW) effect because of the different 
curvature and cosmological constant content of the models. 
From the practical point of view, however, it is still very 
difficult to break the degeneracy, since measurements are 
limited by "cosmic variance" on those scales, and because 
of the possible contribution of gravitational waves. 

Allowing AiV to vary, but keeping constant the other 
3 parameters Ub, z eq , and 1Z, we obtain nearly degenerate 
power spectra which we plot in Fig. 1, normalized to the 
first acoustic peak. The degeneracy in the acoustic peaks 
region is now slightly spoiled by the variation of the ratio 
Q,~ l /Q. re i: the different radiation content at decoupling in- 
duces a larger (for AN > 0) early ISW effect, which boosts 
the height of the first peak with respect to the other acoustic 
peaks. Nevertheless, it is still impossible to distinguish be- 
tween the different models with present CMB measurements 
and without external priors. Furthermore, a slight change in 
the scalar spectral index, n s , can reproduce a perfect degen- 
eracy up to the third peak. 

The main result of this is that, even with a measure- 
ment of the first 3 peaks in the angular spectrum, it is im- 
possible to put bounds on u) re i alone, even when fixing other 
parameters such as u)b- Furthermore, since the degeneracy 
is mainly in z eq , the constraints on ui m from CMB are also 
affected (see section ^J). 

In Fig. 2 we plot the shift parameter 1Z as a function of 
AiV, while fixing fi m = 0.3 and S1a = 0.7. Increasing AiV 
moves the peaks to smaller angular scales, even though the 
dependence of the shift parameter on AiV is rather mild. 
In order to compensate this effect, one has to change the 
curvature by increasing Q m and Qa- We therefore conclude 
that the present bounds on the curvature of the universe are 
weakly affected by AiV. Nevertheless, when a positive (neg- 
ative) AiV is included in the analysis, the preferred models 
are shifted toward closed (open) universes. 



4 CMB ANALYSIS 

In this section, we compare the recent CMB observations 
with a set of models with cosmological parameters sampled 
as follows: 0.1 < Cl m < 1.0, 0.1 < fi re i/fi re ! (AiV = 0) < 3, 



0.015 <Q b < 0.2; < n A < TO and 0.40 <h< 0.95. We 
vary the spectral index of the primordial density perturba- 
tions within the range n s — 0.50, 1.50 and we re-scale 
the fluctuation amplitude by a pre-factor Cio, in units of 
C? OBE . We also restrict our analysis to purely adiabatic, 
flat models (Qtot = 1) and we add an external Gaussian 
prior on the Hubble parameter h — 0.65 ± 0.2. 

The theoretical models are computed using the pub- 
licly available CMBFAST program (Seljak & Zaldarriaga 1996) 
and are compared with the recent BOOMERanG-98, DASI 
and MAXIMA- 1 results. The power spectra from these ex- 
periments were estimated in 19, 9 and 13 bins respec- 
tively, spanning the range 25 < I < 1100. We approxi- 
mate the experimental signal C B X inside the bin to be a 
Gaussian variable, and we compute the corresponding theo- 
retical value Cg 1 by convolving the spectra computed by 
CMBFAST with the respective window functions. When 
the window functions are not available, as in the case of 
Boomerang-98, we use top-hat window functions. The like- 
lihood for a given cosmological model is then defined by 
-21nL = (Cf - C B x )M BB ,(Cf, - C&) where Cf (Cj?) 
is the theoretical (experimental) band power and M BB i is 
the Gaussian curvature of the likelihood matrix at the peak. 
We consider 10%, 4% and 4% Gaussian distributed calibra- 
tion errors (in fi K) for the BOOMERanG-98, DASI and 
MAXIMA- 1 experiments respectively. We also include the 
COBE data using Lloyd Knox's RADPack packages. 

In order to show the effect of the inclusion of uj re i on 
the estimation of the other parameters, we plot likelihood 
contours in the ui re i — iv m , oj re i — tub, w re ; — n a planes. 

Proceeding as in Melchiorri et al. (2000), we calculate a 
likelihood contour in those planes by finding the remaining 
'nuisance' parameters that maximize it. We then define our 
68%, 95% and 99% confidence levels to be where the likeli- 
hood falls to 0.32, 0.05 and 0.01 of its peak value, as would 
be the case for a 2-dimensional Gaussian. 

In Fig. 3 we plot the likelihood contours for cu re i vs 
u) m ,u)b and n s (top to bottom). As we can see, ui re i is very 
weakly constrained to be in the range 1 < ui re i/u} re i(AN = 
0) < 1.9 at 1— a in all the plots. The degeneracy between u re i 
and u>m is evident in the top panel of Fig. 3. Increasing ui re i 
shifts the epoch of equality and this can be compensated 
only by a corresponding increase in u) m . It is interesting 
to note that even if we are restricting our analysis to flat 
models, the degeneracy is still there and that the bounds 
on u m are strongly affected. We find ui m — 0.2 ± 0.1, to 
be compared with u m = 0.13 ± 0.04 when AN is kept to 
zero. It is important to realize that these bounds on u re i 
appear because of our prior on h and because we consider 
flat models. When one allows h as a free parameter and 
any value for O m , then the degeneracy is almost complete 
and there are no bounds on cu re i . In the middle and bottom 
panel of Fig. 3 we plot the likelihood contours for uib and 
n s . As we can see, these parameters are not strongly affected 
by the inclusion of Lu re i. The bound on u>b, in particular, is 
completely unaffected by uj re i. There is however, a small 
correlation between u re i and n s : the boost of the first peak 
induced by the ISW effect can be compensated (at least up 
to the third peak) by a small change in n s . 

Since the degeneracy is mainly in z eq , it is useful to esti- 
mate the constraints we can put on this variable. In Fig. 4 we 
plot the likelihood contours on z eq by using the marginal- 
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ization/maximization algorithm described above. By inte- 
gration of this probability distribution function we obtain 
z eq — 31001|qq at 68% c.L, i.e. a late-time equality, in agree- 



ment with a low-density universe. 

External constraints. It is interesting to investigate how 
well constraints from independent non-CMB datasets can 
break the above degeneracy between tj re ; and u m . The su- 
pernovae luminosity distance is weakly dependent on u) re i 
(see however Zentner & Walker 2001), and the bounds ob- 
tained on Q. m can be used to break the CMB degeneracy. 
Including the SN-Ia constraints on the fi m — Qa plane, 
0.8O m - 0.6fi A = -0.2 ± 0.1 (Perlmutter et al. 1999), we 
find u> re i/uj re t(AN = 0) = 1.12 ( L'o 5 42 at the 2 — a confidence 
level. 

It is also worthwile to include constraints from galaxy 
clustering and local cluster abundances. The shape of the 
matter power spectrum in the linear regime for galaxy clus- 
tering can be characterized by the shape parameter V ~ 
Qmh/y/(l + 0.135AiV)e- (ni > (1+v ^ r/n '" ) - ' 06) . From the ob- 
served data one has roughly (see e.g., Bond & Jafle 1998) 
0.15 < T + (n, - l)/2 < 0.3. 

The degeneracy between u) m and ui re i in the CMB can- 
not be broken trivially by inclusion of large-scale structure 
data, because a similar degeneracy affects the LSS data as 
well (see e.g. Hu et al 1999). However, the geometrical de- 
generacy is lifted in the matter power spectrum, and accu- 
rate measurements of galaxy clustering at very large scales 
can distinguish between various models. This is exempli- 
fied in the bottom panel of Fig. 1, where we plot 3 matter 
power spectra with the same cosmological parameters as in 
the top panel, togheter with the decorrelated matter power 
spectrum obtained from the PSCz survey. 

The inclusion of the above (conservative) value on T 
gives uj re i/uj re i(AN = 0) = 1.40_o756j that is less restrictive 
than the one obtained with the SN-Ia prior. 

A better constraint can be obtained by including a prior 
on the variance of matter perturbations over a sphere of 
size 8/i -1 Mpc, derived from cluster abundance observa- 
tions. Comparing with ag = (0.55 ± 0.05)f2m°' 47 , we obtain 
Urel J Wei(AjV = 0) = 1.27°'o 5 4 3 , again at the 2 — a confi- 
dence level. 
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5 FORECAST FOR MAP AND PLANCK 

In this section we perform a Fisher matrix analysis in order 
to estimate the precision with which forthcoming satellite 
experiments will be able to constrain the parameter z £q . 

Fisher matrix. Using C(s) to denote the likelihood func- 
tion for the parameter set s and expanding In C to quadratic 
order about the maximum defined by the reference model 
parameters so, one obtains 

C w £(s ) exp(— i FijSsiSsj) 



where the Fisher matrix Fij is given by the expression 



E 



dC e dC e 



(ACf) 2 dsi d Sj 



(6) 



and £ m ax is the maximum multipole number accessible to 
the experiment. The quantity ACi is the standard devia- 
tion on the estimate of Ci, which takes into account both 
cosmic variance and the expected error of the experimental 
apparatus and is given by 

(AC,) 2 » fna t \ t (C e +B; 2 ) 2 , (7) 



B'r 



E 



(21 + l)f sky 



(8) 



(Knox 1995; Efstathiou & Bond 1999), for an experiment 
with N channels (denoted by a subscript c), angular resolu- 
tion (FWHM) 9 C , sensitivity a c per resolution element and 
with a sky coverage f s ky The inverse weight per solid an- 
gle is Wc 1 = (o~ c 9c)~ 2 and £ c = \/8 In 2/8 c is the width of 
the beam, assuming a Gaussian profile. If the initial fluc- 
tuations are Gaussian and a uniform prior is assumed, one 
finds that the covariance mat rix is given by th e inverse of 
the Fisher matrix, C = F" 1 ( fcond et al. 1997| ). The stan- 
dard deviation for the parameter s, (with marginalization 
over all other parameters) is therefore given by Oi — \JCu. 
This approximation is rigorously valid only in the vicinity of 
the maximum of the likelihood function, but it has proved to 
give useful insight even for large values of s — so (Efstathiou 
& Bond 1999; Efstathiou 2001). The main advantage of the 
Fisher matrix approach when compared to an exact like- 
lihood analysis is that for m cosmological parameters the 
former requires only the evaluation of m + 1 power spectra. 
Therefore the computational effort is vanishingly small with 
respect to the one necessary for a full likelihood analysis of 
the parameter space. 

Table 1 summarizes the experimental parameters for 
MAP and Planck we have used in the analysis. For both 
experiment we have taken f s k y = 0.50. These values are 
indicative of the expected performance of the experimental 
apparatus, but the actual values may be somewhat different, 
especially for the Planck satellite. 

Cosmological parameters. The validity of the Fisher ma- 
trix analysis depends on the chosen parameter set, as well 
as on the point so at which the likelihood function is sup- 
posed to reach its maximum. We use the following 9 di- 
mensional parameter set: LUb, u) c , loa, TZ, z eq ,n s ,n t ,r, Q. Here 
n s , m are the scalar and tensor spectral indices respectively, 
r = Cj/Cf is the tensor to scalar ratio at the quadrupole, 
and Q =< 1(1 + l)Ct > 1 ' /2 denotes the overall normaliza- 
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Fig. 5 — Derivatives of with respect to the 9 parameters eval- 
uated at the reference model described in the text. The derivative 
dC'i/duiA has been set to for £ > 200 in order to suppress the ef- 
fect of numerical errors, thus taking into account the geometrical 
degeneracy. Fig. (a): 0.1 -dCi/duJt, (red, solid); dCi/dz eq (green, 
dotted); dC e /duJA (blue, dashed). Fig. (b): 10 • dC e /dr (yellow, 
solid); dCi /du> c (green, dotted); dCi/dlZ (black, dashed). Fig. 
(c): dC(/dn s (light blue, solid); 10 ■ dCi/dnt (magenta, dotted); 
dC'i/dQ (black, dashed). 
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0.032 
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3.170 
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5n 3 


0.149 


0.013 
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8r/r 


5.222 
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SQ 


0.012 
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Table 2 — Fisher matrix analysis results: expected 1 — a errors 
for the MAP and Planck satellites. See the text for details and 
discussion. 

racy, while Planck will pinpoint it down to 2% or so. From 
u) re i = (uib + u)c)/z eq it follows that the energy density of 
relativistic particles, ui re i, will be determined within 43% by 
MAP and 3% by Planck. This translates into an impossibil- 
ity for MAP of measuring the effective number of relativistic 
species (AN e ff ~ 3.17 at la), while Planck will be able to 
track it down to AN e jj ~ 0.24. As for the other parame- 
ters, while the acoustic peaks' position (through the value 
of TV) and the matter content of the universe can be de- 
termined by Planck with high accuracy (of the order of or 
less than one percent), the cosmological constant remains 
(with CMB data only) almost undetermined, because of the 
effect of the geometrical degeneracy. The scalar spectral in- 
dex n s and the overall normalization will be well constrained 
already by MAP (within 15% and 1%, respectively), while 
because of the reasons explained above the tensor spectral 
index n t and the tensor contribution r will remain largely 
unconstrained by both experiments. Generally, an improve- 
ment of a factor 10 is to be expected between MAP and 
Planck in the determination of most cosmological parame- 
ters. Our analysis considers temperature information only, 
while it is well known that inclusion of polarization measure- 
ments greatly improves determination of the tensor mode 
parameters (see e.g. Bucher et al. 2000). For Planck, which 
will have polarization measurement capabilities, this will be 
of great importance. 

A Fisher matrix analysis for AN e f f was previously per- 
formed by Lopez et al. (1999) and repeated by Kinney & Ri- 
otto (1999) (with the equivalent chemical potential £), and a 
strong degeneracy was found between N e f j , h and Qa , and 
to lesser extent with Q b . We have seen here that the degen- 
eracy really is between u} re i,uj m and n, and the degeneracy 
previously observed is thus explained because they consid- 
ered flat models, where a change in Oa is equivalent to a 
change in uj m , u) m = (1 — ^a — Q. b )h 2 . The results of this 
paper, on how precisely the future satellite missions can ex- 
tract the relativistic energy density, can be translated into 
approximately AN eff = 3.17 (£ = 2.4) and AN eff = 0.24 
(£ = 0.73) for MAP and Planck respectively. 



Table 1 — Experimental parameters used in the Fisher matrix 
analysis. 

tion, where the mean is taken over the multipole range 
accessible to the experiment. We choose to use the shift 
parameter 1Z because this takes into a ccount the geomet 



rical degeneracy between D.a and Q,k (Efstathiou & Bond 
ur purely adiabatic reference model has parame- 
ters: oj b = 0.0200 (fi 6 = 0.0473), lj c = 0.1067 (fi c = 0.2527), 
'jja = 0.2957 (fl A = 0.7000), (h = 0.65), 1Z = 0.953, 
Ze q = 3045, n a = 1.00, n t = 0.00 , r = 0.10, Q = 1.00. 
This is a fiducial, concordance model, which we believe is in 
good agreement with most recent determinations of the cos- 
mological parameters (flat universe, scale invariant spectral 
index, BBN compatible baryon content, large cosmological 
constant). Furthermore, we allow for a modest, 10% tensor 
contribution at the quadrupole in order to be able to include 
tensor modes in the Fisher matrix analysis. 

We plot the derivatives of Ce with respect to the differ- 
ent parameters in Fig. 5. Generally, we remark that deriva- 
tives with respect to the combination of parameters describ- 
ing the matter content of the universe (u> b and U) c , TZ, z eq ) 
are large in the acoustic peaks region, i > 100, while deriva- 
tives with respect to parameters describing the tensor con- 
tribution (nt, r) are important in the large angular scale 
region. Since measurements in this region are cosmic vari- 
ance limited, we expect uncertainties in the latter set of 
parameters to be large regardless of the details of the ex- 
periment. The curve for dCe/dQ is of course identical to 
the Ct's themselves. The cosmological constant is a notable 
exception: variation in the value of uja keeping all other pa- 
rameters fixed produces a perfect degeneracy in the acoustic 
peaks region. Therefore we expect the derivative dCe/duiA 
to be in this region. Small numerical errors in the compu- 
tation of the spectra, however, artificially spoil this degener- 
acy, erroneously leading to smaller predicted uncertainties. 
In order to suppress this effect, we set dCe/duuA = for 
I > 200. From eq. (^) we see that a large absolute value 
of dCi/dsi leads to a large Fa and therefore to a smaller 
1 — a error (roughly neglecting non-diagonal contributions). 
If the derivative along s% can be approximated as a linear 
combination of the others, however, then the correspond- 
ing directions in parameter space will be degenerate, and 
the expected error will be important. This is the case for 
mild, featureless derivatives as dCe/dr, while wild changing 
derivatives (such as dCi / dTZ) induce smaller errors in the de- 
termination of the corresponding parameter. Therefore the 
choice of the parameter set is very important in order to 
correctly predict the standard errors of the experiment. 

Error forecast. Table 2 shows the results of our analy- 
sis for the expected 1 — a error. Determination of the red- 
shift of equality can be achieved by MAP with 23% accu- 
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6 CONCLUSIONS 

In this paper, we have examined the effect of varying the 
background of relativistic particles on the cosmological pa- 
rameters derived from CMB observations. We have found 
that the present constraints on the overall curvature, flk, 
and tilt of primordial fluctuations, n s , are slightly affected 
by the inclusion of this background. However, we have found 
a relevant degeneracy with the amount of non relativistic 
matter ui m . Even with relatively strong external priors (flat- 
ness, h = 0.65±0.2, age > 10 Gyrs) the present CMB bound 
(95% c.l.) 0.1 < uj m < 0.2 spreads to 0.05 < w m < 0.45 when 
variations in uj re i are allowed. Specifically, without priors on 
uim (through flatness, h, etc) no bounds on N e ff can be 
obtained. 

Another fundamental point bears on the identification 
of the best choice of parameters, i. e. parameter combina- 
tions which can be unambigously extracted from CMB data. 
It is of the greatest importance to realize which parameter 
set is least plagued by degeneracy problems, i. e. which di- 
rections in parameter space are non- flat. In the well known 
case of the geometrical degeneracy, the shift parameter 1Z 
can be determined with very high precision by measuring 
the position of the peaks. The curvature and the Hubble 
parameters, however, are almost flat directions in parame- 
ter space, and therefore are not ideal variables for extraction 
from CMB data. In this work, we have pointed out that an 
analogous situation exists for z eq , ui re i and ui m - In fact, z eq is 
well determined because it measures the physical distance to 
equality time, while on the contrary uj m is a rather ill-suited 
variable for CMB data, since it suffers from degeneracy with 
u>rei (at least up to the third acoustic peak). 

Fortunately, as we saw in the last section, the matter 
- radiation degeneracy in the CMB data is present only up 
to the third peak and future space missions like Planck will 
be able to determine separately the amount of matter and 
radiation in the universe. 
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